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Abstract: 

Recently a non-supersymmetric analog of Veneziano-Yankielowicz (VY) effective La- 
grangian has been proposed and applied for the analysis of the 9 dependence in pure 
Yang-Mills theory. This effective Lagrangian is similar in many respects to the VY con- 
struction and, in particular, exhibits a kind of low energy holomorphy which is absent 
in the full YM theory. Here we incorporate a heavy fermion into this effective theory by 
using the "integrating in" technique. We find that, in terms of this extended theory, holo- 
morphy of the effective Lagrangian for pure YM theory naturally implies a holomorphic 
dependence on the heavy fermion mass. It is shown that this analysis fixes, under certain 
assumptions, a dimensionless parameter which enters the effective Lagrangian and deter- 
mines the number of nondegenerate vacuum sectors in pure YM theory. We also compare 
our results for the vacuum structure and 9 dependence to those obtained recently by 
Witten on the basis of AdS/CFT correspondence. 



1 Introduction 



The remarkable progress made over the past years in understanding the nonperturbative 
properties of supersymmetric (SUSY) N=l and N=2 theories owes a lot to holomorphy. 
Supersymmetry requires a superpotential W e ff(X i ,g i ) of the effective low energy theory 
to be a holomorphic function of the light fields and the coupling constants <?j, and thus 
powerfully constraints the large distance dynamics of SUSY theories (see e.g. for 
a review). The famous Veneziano-Yankielowicz (VY) effective Lagrangian ||, obtained 
from the anomalous Ward identities of SUSY YM theory, is holomorphic in its fields and 
parameters[|. Moreover, the effective VY potentials for theories with different numbers of 
matter fields (e.g. ones with Nf and Nf — 1 fermion flavors) are related to each other by 
holomorphic decoupling relations More recently, an operation inverse to integrating 
out of a heavy fermion was suggested under the name of "integrating in" procedure Q. 
The integrating in technique provides a powerful method to obtain relations between 
parameters and functional forms of superpotentials for theories with different matter 
contents. In particular, it was shown in Q that Affleck-Dine-Seiberg superpotential |J 
can be obtained from VY Lagrangian for SUSY gluodynamics by the integrating in/out 
procedure. The integration in technique was used to study the phase structure of SUSY 
theories 0, and to calculate numerical parameters in different models J7|. 

Recently a non-supersymmetric analog of the VY effective Lagrangian has been pro- 
posed H using an infinite series of anomalous Ward identities for YM theory. This con- 
struction has some striking similarities to its supersymmetric counterpart. Specifically, 
we have found that this effective Lagrangian (more precisely, effective potential) possesses 
both a "dynamical" part, which is similar to the original VY form ||, and a "topological" 
part, which is analogous to an improvement of the VY Lagrangian, suggested recently 
0. The picture of the physical 9 dependence in pure YM theory, following from the 
analysis of this effective Lagrangian, is rather similar to the one found for SUSY YM 
theory |10|, ||: the correct 2n periodicity in 6 is recovered when a set of disconnected 
vacua is taken into account. The difference from the supersymmetric case is in the ab- 
sence of degeneracy between the vacua. As a result, all vacua but one with a lowest (for 
a fixed 9) energy contribute zero to the partition function in the thermodynamic limit 
||. The number of vacua can be found from explicit dynamical calculations only. For 
the SUSY case, such calculations can be done [|K|, while no similar technique is known 
in the non-supersymmetric case. The information on the number of vacua was coded in 
the effective Lagrangian of Ref. |§ in a dimensionless number £ ~ N~ l (see Eq.(^) below). 
(On general grounds, the number of vacua in YM theory should be proportional to the 
number of colors N c .) 

Another remarkable property of the effective potential of Ref. is its holomorphic 
structure. At first sight, a claim of holomorphy of an effective Lagrangian for YM theory 
may sound suspicious, as the full theory is not seen to possess any holomorphy, in contrast 
to supersymmetric models. As will be discussed in detail below, there is no contradiction 
here: holomorphy is the property of a nonperturbative effective potential (defined within 
a particular regularization scheme) which is fixed in this approach by the anomalous 

x It is well known that the VY effective Lagrangian is not a genuine Wilsonian effective Lagrangian for 
light degrees of freedom, but rather has a different meaning, see Sect. 2. 
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Ward identities and some additional arguments, and describes the large distance physics 
only. We do not expect that kinetic terms and/or perturbative contributions would be 
holomorphic as well, but they are irrelevant anyway for our purposes. 

In this paper we analyse holomorphy of the effective Lagrangian of Ref. || using the 
integration in technique [|J. By the integrating in procedure we obtain an effective La- 
grangian for the YM field interacting with a heavy fermion. It will be shown that holo- 
morphy of the effective Lagrangian for pure YM theory corresponds to a holomorphic 
dependence on a fermion mass in an extended theory including the heavy fermion. When 
gluodynamics is defined as a low energy limit of the latter theory, the integrating in 
procedure fixes (under certain plausible assumptions specified below) the aforementioned 
parameter £, and thus allows one to find the number of different nondegenerate 9 vacuum 
sectors H for pure YM theory. As will be shown below, this approach results in the 
same value £ = 4/(36) which was obtained in [] L| with a different method which implies, 
however, the same regularization prescription. 

Our presentation is organized as follows. In Sect. 2 we recall the construction || of the 
effective Lagrangian for pure gluodynamics. In Sect. 3 a new effective Lagrangian for the 
theory including a heavy fermion is obtained by the integration in technique, assuming 
the standard form of the fermion mass term and preservation of holomorphy under the 
procedure of integrating in/out. We show that the parameter £ can be fixed by comparing 
the holomorphic properties of two effective Lagrangians. The matching of two effective 
theories is further discussed in Sect. 4 from the point of view of their global, "topological" , 
properties. We also study in this section a connection of the present approach with the 



previous analysis of Ref. ||1 1|| . In Sect. 5 we discuss our results and compare them with 
different approaches to the problems of interest. In particular, we argue that our picture 
of the vacuum structure and 9 dependence is in qualitative agreement with that recently 
obtained by Witten |T3| in the limit N c — > oo within a different approach. The final Sect. 6 
contains our conclusions. 



2 Effective Lagrangian for gluodynamics 

The purpose of this section is to describe an effective Lagrangian for YM theory, which 
was constructed in our previous paper ||. Before proceeding with the presentation, we 
would like to pause for a comment on the meaning of this effective Lagrangian. As there 
exist no Goldstone bosons in pure YM theory, no Wilsonian effective Lagrangian, which 
would correspond to integrating out heavy modes, can be constructed for gluodynamics. 
Instead, one speaks in this case of an effective Lagrangian as a generating functional for 
vertex functions of the composite fields G 2 and GG. Moreover, only the potential part 
of this Lagrangian can be found as it corresponds to zero momentum n-point functions 
of G 2 , GG, fixed by arguments appealing to renormalizability of the theory. The kinetic 
part is not fixed in this way. Thus, such an effective Lagrangian is not very useful for 
calculating the S-matrix, but is perfectly suitable for addressing the vacuum properties^. 
Specifically, space-time independent fields are amenable to a study within this framework. 

2 For SUSY theories, Veneziano-Yankielowicz effective Lagrangian H has the same meaning, see M. 
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Analogously to the case of SUSY theories, an effective Lagrangian for pure YM the- 
ory is constructed using an infinite series of anomalous Ward identities which serve as 
matching conditions ensuring consistency of the large distance properties of the theory 
with its small distance behavior fixed by renormalizability and asymptotic freedom. The 
complete set of two-point correlation functions is 

lim i J dxe^(0\T{^G 2 (x) ^G 2 (0)}|0> = "4<^G 2 > , (1) 

lim i J dxe^(0\T j^W) ^GG(O) j |0> = ~M^GG) , (2) 

lim i / dxe^(0\T{^GG(x)^GG(0)\\0)=e(^G*) . (3) 
9^0 J {Sir 87T J Aa s 

Here £ is a dimensionless parameter which will be assumed to be a rational number. (An 
irrational value of £ would presumably produce a non-differentiable 9 dependence for YM 
theory.) Multi-point correlation functions of the fields G 2 ^ and G^ V G^ U are obtained by 
differentiating Eqs.([|-|3|) with respect to 9 and l/g 2 , see below. In what follows we use 
the one-loop /3-function /3(a s ) = — 6a 2 /(27r) where b = (ll/3)iV c , though most of the 
discussion below can also be formulated with formally keeping the full /3-function. 

A few comments on these Ward identities are in ordeiQ. All correlation functions 
(0),(0),(0) are defined via Wick type of the T-product, i.e. by the differentiation of the 
path integral in respect to corresponding parameters, see below. Perturbative contribu- 
tions to the conformal anomaly matrix element (—ba s /(8Tr)G 2 ) in Eqs. ([I]) and ([|) are 
subtracted to any finite order in a s by definition. In this case, its dependence on the bare 
coupling constant g corresponding to the cut-off scale Mr is fixed by the dimensional 
transmutation formula 

4 

= °i A yM • (4) 

Analogously, (a s /(8ir)GG) = c 2 Ay M . Here numerical constants Ci,c 2 are independent 
of go, but depend on the vacuum angle 9 which is allowed to be non-zero in Eqs.(|l]-^D. 
We note that the constants c\,c-i depend on a particular regularization scheme used to 
define the nonperturbative vacuum condensates. However, once specified, the VEV (§) 
determines all zero momentum correlation functions of f3(a s )/(4:a s )G 2 , with perturbative 
tails subtracted. The Ward identities ([[]), (|2|) then follow by the differentiation of the 
above expressions with respect to 1/g 2 . They were derived long ago by Novikov, Shifman, 
Vainshtein and Zakharov (NSVZ) [0. By derivation, the two-point functions ([lfj^) do 
not contain perturbative contributions^. As was discussed in detail in Ref.||], the zero 

3 Although we will occasionally call Eqs. (0)-(||) the Ward identities, it should be mentioned that Eq. 
is not precisely a Ward identity, but rather is a scheme dependent relation involving an unspecified 

at this stage parameter £. 

4 The same result (|l|) was obtained in [fl3[ using canonical methods with Pauli-Villars regularization. 
To one loop order in regulator fields, it was found that perturbative contributions add an identity to both 
sides of Eq. (0) . The absence of perturbative contributions to the correlation functions (||) , (||) is obvious. 
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momentum correlation functions ([T]-|3|) and their n-point generalizations are generated by 
the differentiation of \og(Z/ZpT) with respect to 1/g 2 and 9, where 

Z{9) = Z PT exp {-iVE v (9)} = Z PT exp |-iV<0| - ^C 2 |0) e | . (5) 

Here the perturbatively defined partition function Zp? does not depend on 9 and absorbs 
perturbative contributions to the conformal anomaly. The dimensionless parameter £ in 
Eq.(^|) is related to the 9 dependence of the vacuum energy E v {6) in Eq.(|J): 

E v {9) = E v (0) f{9) , f(9) = 1 - 2£ 2 9 2 + ■ ■ ■ . (6) 

Let us introduce complex linear combinations of the composite fields 

In terms of these combinations, the Ward identities (0-0) take particularly simple forms 
(for an arbitrary value of the vacuum angle 9): 

lim i f dxe iqx (0\T{H(x) #(0)}|0) = -4(H) , 

q^O J 

limi J dxe tqx (0\T{H{x) #(0)}|0) = -4(F), (8) 

lim i f dxe iqx (Q\T{H(x) H(0)}\0) = . 
q^o J 

It can be seen that the n-point zero momentum correlation function of the operator H 
equals (— 4) n ~ l (H). Multi-point correlation functions of the operator H are analogously 
expressed in terms of its vacuum expectation value (H). At the same time, it is easy to 
check that the decoupling of the fields H and H holds for arbitrary n-point functions of 
if, H. This is the origin of holomorphy of an effective Lagrangian for YM theoryQ. 

An effective low energy Lagrangian (more precisely, effective potential) is now con- 
structed as the (Legendre transform of) generating functional for zero momentum cor- 
relation functions of the marginal operators G^G^ and G^G^, such as Eq.(H) and 
their n-point generalizations. It is a function of effective zero momentum fields h, h which 
describe the vacuum expectation values of the fields H, H : 

J dxh = (J dx H) , J dxh = (J dx H) . (9) 



5 To avoid possible misunderstanding, we note that the right hand side of the last equation in Q does 
contain perturbative contributions proportional to regular powers of a s . However, they are irrelevant 
for our purposes. The decoupling of the fields H and H holds at the level of nonperturbative 0(e -1 /" s ) 
effects. Holomorphy of an effective potential for YM theory has the same status. Thus, in contrast to the 
supersymmctric case where holomorphy is an exact property of the effective supcrpotcntial, in the present 
case it only refers to a "nonperturbative" effective potential which does not include perturbative effects 
to any finite order in a s . We assume that perturbative and nonperturbative effects can be separated, at 
least in principle or/and by some suitable convention [ fj"l| ||. 
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Omitting details of the derivation, which can be found in ||, we give the final answer for 
the improved effective potential F(h, h): 

e ~ iVFm = %3^{~^{ hLo9 ^ +jiL ° 9 i) 

+ nrV ik+- — — — , 10 

\ P 2lX J 2l J 

where the constants C, C can be set real and expressed in terms of the vacuum energy 
at 6 = 0, C = C = —2eE v (0), see Eq.fl^), and V is the 4- volume. The holomorphic 
structure of the effective potential is explicit in Eq.(|H]). The integer numbers p and q are 
relatively prime. As was shown in Ref. 0, they are related to the parameter £ introduced 
in Eq.(|3]), 2£ = q/p. The symbol Log in Eq.fllOD stands for the principal branch of the 
logarithm. The effective potential flTOl) produces an infinite series of anomalous WI's. By 
construction, it is a periodic function of the vacuum angle 9. The effective potential (|T0|) 
is suitable for a study of the YM vacuum as described above. 

The double sum over the integers n, k in Eq. (|3|) appears as a resolution of an ambiguity 
of the effective potential as defined from the anomalous WI's. As was discussed in f|, 
this ambiguity is due to the fact that any particular branch of the mult i- valued function 

log z v l q = - Log z + 2m(n + k~) , n = 0,±l,... ; k = 0, 1, . . . , q - 1 (11) 
q q 

corresponding to some fixed values of n, k, satisfies the anomalous WI's. However, without 
the summation over the integers n, k in Eq . (|i0|) , the effective potential would be multi- 
valued and unbounded from below. An analogous problem arises with the original VY 
effective Lagrangian. It was cured by Kovner and Shifman in by a similar prescription 
of summation over all branches of the multi-valued VY superpotential. Moreover, the 
whole structure of Eq. ([H]) is rather similar to that of the (amended) VY effective potential. 
Namely, it contains both the "dynamical" and "topological" parts (the first and the second 
terms in the exponent, respectively). The "dynamical" part of the effective potential fll0|) 
is similar to the VY || potential ~ Slog(S/A) Nc (here S is an anomaly superfield), 
while the "topological" part is akin to the improvement of the VY effective potential. 
Similarly to the supersymmetric case, the infinite sum over n reflects the summation 
over all integer topological charges in the original YM theory. The difference of our 
case from that of supersymmetric YM theory is that an effective potential of the form 
(l/A^)01og(0/A) Ar , as in the SUSY case, implies a simpler form of the "topological" term 
~ 2nin/N (cf) — 0) with only one "topological number" n which specifies the particular 
branch of the multi-valued logarithm. In our case, we allow for a more general situation 
when the parameter N is a rational number N = p/q. In this case we have two integer 
valued "topological numbers" n and k, specifying the branches of the logarithm and 
rational function, respectively. Our choice is related to the fact that some proposals to fix 
the values of p, q suggest that q ^ 1. One may expect the integers p and q are related to a 
discrete symmetry surviving the anomaly, which may not be directly visible in the original 
fundamental Lagrangian. Our purpose in this paper is to try to find these numbers by 
analysing the properties of the effective potential ([H]). 
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It should be stressed that the improved effective potential (P]) contains more infor- 
mation in comparison to that present in the anomalous Ward identities just due to the 
presence of the "topological" part in Eq.([H]). Without this term Eq.flIDD would merely be 
a kinematical reformulation of the content of anomalous Ward identities for YM theory. 
The reason is that the latter refer, as usual, to the infinite volume (thermodynamic) limit 
of the theory, where only one state of a lowest energy (for 9 fixed) survives. This state 
corresponds to one particular branch of the multivalued effective potential in Eq. fllOD . 
At the same time, the very fact of multi-valuedness of the effective potential implies that 
there are other vacua which should all be taken into consideration when 9 is varied. When 
summing over the integers n , k, we keep track of all (including excited) vacua of the the- 
ory, and simultaneously solve the problems of multi-valuedness and unboundedness from 
below of the "one-branch theory". The most attractive feature of the proposed struc- 
ture of the effective potential (|T0|) is that the same summation over n, k reproduces the 
topological charge quantization and 2ir periodicity in 9 of the original YM theory. 



3 Integrating in the heavy fermion 

In this section the effective Lagrangian (|1(]) for pure YM theory will be analysed from 
a different standpoint. Considering pure gluodynamics as a low energy limit of a theory 
describing the YM field interacting with a heavy fermion, we now wish to construct an 
effective Lagrangian for the latter theory starting from the effective Lagrangian fllPl) . 
Our purpose here is to try to understand, in this framework, holomorphy of the effective 
potential (|T0D and the constraints imposed by it. As will be argued below and in Sect. 4, 
a relation between the holomorphic and "topological" properties of two Lagrangians is 
non-trivial, and allows one to fix the crucial parameter £ = q/(2p) entering Eq. (TlO|) , under 
two plausible assumptions, see below. 

The task of constructing such an effective Lagrangian for the theory with a heavy 
fermion is achieved by using the "integrating in" technique, developed in the content of 
SUSY theories in Ref. j|] and reviewed by Intriligator and Seiberg in 0. The integrating 
in procedure can be viewed as a method of introducing an auxiliary field into the effective 
Lagrangian for pure YM theory. Using the renormalization group properties of the YM 
effective Lagrangian, the latter is extended to include the auxiliary field U, which will be 
later on identified with the chiral combination l§> r of a heavy fermion. 

To conform with the notation and terminology of Ref. ||], we will call pure YM theory 
and the theory with a heavy fermion the d-theory (from "downstairs" ) and the u-theory 
(from "upstairs"), respectively. The effective potential of the d-theory is then Wd + Wj~ 
with 

w *i h u*Y' = li h ^>* (12) 

Ap \C J Ap (cA Y M) p/q 

(here c is a dimensionless numerical coefficient), and the summation over all branches 
of the logarithm in the partition function is implied. In this section Eq.(^) will be 
understood as representing a branch (section) of the multi- valued effective potential, which 
corresponds to a lowest energy state for small 9 n. As was shown in ||, this section 
corresponds to the principal branch of the rational function in the logarithm in Eq. (|P2|) . 
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In this case the vacuum expectation value (H) depends on 9 as follows: 



(H) e = (H) e 2 W (13) 

We now want to relate jll], || the dimensional transmutation parameter Aym of pure 
YM theory to the scale parameter Aqcd of the u-theory including a heavy quark of mass 
m ^> Aqcd, Aym- We assume both parameters to be defined in the MS scheme, in which 
no threshold factors arise in corresponding matching conditions. The matching condition 
then follows from the standard one-loop relations 

Aym = Moexp (- ^ (M „) ) • = b = ^ JV C , 

A QCD = Moexp (- bQc ^ (Uo) ) , » Q ao = ^N e - 2 -, (14) 

and the requirement that the coupling constants of the d- and u- theories coincide at the 
decoupling scale M = m. We obtain 



/ m 2 \ 4/(36) 

A-YM = ^QCD T2 ' (^) 



As was explained in Refs. |L0|, @], Eq. ([l5|) reflects the fact that, for fixed Aqcd, the scale 



parameter Aym characterizes the low energy theory surviving below the scale m, and thus 
depends on m. In this sense, the constant in the logarithm in Eq . (|12|) also depends on m 



8p/(3bq) 

(c^YMY /q = (cAqcd)^ [ ~T~~ ) • (16) 




Following Ref . [Q , we now wish to consider (a particular branch of) the effective potential 
(ITS]) as the result of integrating out the auxiliary field U in the new effective potential 
W = W u — mU which corresponds to the u-theory: Wd(h, m) = W(h, m, (U)), or 

W d =[W u -mU] {u) , (17) 

where (U) is a solution of the classical equation of motion for the auxiliary field U : 

— -m = 0. (18) 

(At this stage, Eq.(|17|) is merely a definition; later on, the term mU will be identified 
with the fermion mass term in the effective Lagrangian.) Let us note that, according to 
Eq . (|17|) , Wd should depend holomorphically on (U). Our assumption is that this is only 
possible if an effective potential W of the u-theory is itself holomorphic in the field U. In 
our opinion, this assumption appears to be quite reasonable^]. Furthermore, one can see 
that Eqs. ([□]), (0) actually define the potential Wd as the Legendre transform of W u . 

6 We are not aware of any counter-example where a holomorphic d-potential would be obtained from 
a non-holomorphic u-potcntial by the integrating out procedure. 
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Therefore we can find the unknown function W u from the known potential W d by the 
inverse Legendre transform: 

W u = [W d + mU] {m) , (19) 

where (m) solves the equation 



_d_ 

dm 



(W d + mU) = 0. 



(20) 



Eq.(P0|) can be considered as an equation of motion for the auxiliary "field" m. It is 
important to note that Eqs. (|17] - ^0|) imply that m should be treated as a complex 
parameter to preserve the holomorphic structure of Eq . (|l~2|) . When substituted in Eq. ([19]) , 
a solution (m) of Eq. (|20|) defines the potential W u (h, U, (m)). When this function is found, 
the effective potential W of the u-theory is defined by the relation 



W(h, U, m) = W u (h, U, (m)) - mU , 

in accord with Eq. (|PT|) . 

The solution of Eq. (|20"D is easy to find using Eqs. (|i~2"D , flTBD 



(21) 



(m) 



2 h 
3~b U 



(22) 



Thus, Eq.( 



) yields 
W u = - 



1 <7 71 

- - h log 
4 p 



'cA 4 



h 



(><d\ ( _2 

J \3b AqcdU 



8p/(3bq)' 



h 



(23) 



Finally, Eq.(^TJ) results in the effective potential of the u-theory 

(l-8/36)p/ 9 / 3fe jj \Sp/(3bqY 



w 



\1 

4 p 



h log 



h 



, UY QCD , 



2 C ^QCD , 



+ 7T7 h - mU . 
3b 



(24) 



We expect this effective Lagrangian to describe YM theory coupled to a heavy quark, 
corresponding to the field U, such that integrating out U brings us back to the effective 
Lagrangian (|l^) for pure gluodynamics. Indeed the equation of motion for the field U 
stemming from the effective potential (|2"4]) reads 



m(U) 



3b 



h. 



(25) 



Inserting this classical vacuum expectation value (VEV) back to Eq.(|2"4]) (i.e. integrating 
out the field U), we reproduce the effective potential of the d-theory, Eq. (|T2"|) . Note that 
as Eq.fl2"o|) should preserve the iV c counting rule, we obtain (h) ~ N% , b ~ A^ c , (U) ~ A^" c . 
The iV c dependence of the VEV (U) is consistent with the identification (U) ~ (^l^r) 
which will be suggested below. 

Analogously to the effective potential (0) of the d-theory, the new potential ( p4"|) is not 
a single- valued function. The single- valuedness should be imposed, as was done in Eq . (|10D , 
by the summation over all branches of W in the partition function. This procedure will be 
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considered in detail in the next section, while here we would like to identify the field U of 
the effective theory with a corresponding operator of the fundamental theory. As is seen 
from (|2"o]), U has dimension 3, and thus should describe the VEV of an operator bilinear 
in the heavy quark fields. Furthermore, as long as m is effectively considered as a complex 
parameter, this operator can only be ^l^r or ^r^l, in accord with the structure of the 
mass term in the underlying fundamental theory. (Here the second assumption of the 
present approach is implicit: we assume the standard form of the fermion mass term in 
the effective potential (53).) To find the exact correspondence, we note that, when the 



VEV of the field h is chosen to correspond to a lowest energy state for 9 < it, Eq . (|25|) 
implies 

(mil + mU) e = -(^ L G 2 ) cos (2£9) , 

(mU - mU) e = -t(-^G 2 ) sin (2f 0) , (26) 
which can be compared with the relations[] between the VEV's in the underlying theory: 
(m* L * B + m* fl ¥ L ) fl = -(^G 2 ) cos (2f 0) + O ( \ 

1Z7T \m z 

(m* L * R - my R q L ) e = -i^ £ (^G 2 ) sin (2£0) + O (-^) . (27) 
Comparing Eqs. ( |2T)| ) and ( P7| ) and using the relation 2^ = g/p ||, we conclude that 

i = 4 > - = r h ; M = <^**> • ( 28 ) 

So p So 

We thus see that the introduction of the heavy quark into the effective theory fixes the 
parameter £ which enters the effective Lagrangian flT2"| ) for pure YM theory. The value 
obtained coincides with the one suggested by us previously within a different method fll 



based on a different assumption. (The correspondence between the two approaches to 
calculate the crucial parameter £ will be further discussed in the next section.) Moreover, 



the comparison of Eqs. (|26|) and (|27 ) shows that the chiral field U corresponds to the 



chiral fermion bilinear of the fundamental Lagrangian: 

mUom^ L ^ R , mU <$7n&R$ L . (29) 

This correspondence between the operators of the effective and underlying theories has 
the same meaning as Eqs.@, i.e. the classical field U describes the VEV of the chiral 
combination ^L^i? of the full theory. 

It may be instructive to discuss the result (E5I) in a slightly different and more intuitive 
way. Using Eqs. ([27]), we can write 



(mU) . m = _i (, + |f) ( ^ GV - - i (i - 



e~ 2 ^, (30) 



7 The equations ( |27| ) follow from the operator product expansions (m\I n I r ) = —(a s / (12n)G 2 ) +0(1 / 'm 2 ) 
(mf i75*) = (a s /(87r)GG) + 0(l/m 2 ), and Eq.(fl|). 
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which is a superposition of the holomorphic (~ exp(2z£#)) and anti-holomorphic (~ 
exp(— 2i£9)) functions. At the same time, the equation of motion (|25|) requires (U) to be 
a holomorphic function, as (h) is such a function. This is only possible, as is seen from 
Eg. (|30"D , when 1 — (36/4)£ = which is equivalent to (|28|). In other words, the requirement 
of holomorphy for the u-theory singles out parameters of the d-theory. 

We would like to pause here to discuss the following issue. The main result of our 
calculations, eq. (p8|) is the direct consequence of two fundamental principles. The first 
principle is a modified definition of the path integral when the prescription of summation 
over all topological sectors is introduced both at the fundamental and effective Lagrangian 
level, see Eq. (|i~0|) . As was discussed earlier, this definition does not change any local 
properties of the theory (in particular, it does not alter the WI's) but drastically changes 
the global properties of the theory. The prescription of summation over the topological 
classes will be further discussed in Sects. 4 and 5. The second fundamental principle we 
adopted from the SUSY theories is the property of holomorphy in the form of a natural 
requirement to have a holomorphic u-theory if a d-theory satisfies this property (and vice 
versa). However, in the SUSY case these two principles are argued to lead to a suspicious 
new chiral invariant vacuum, whose analog certainly can not exist in QCD as we know from 
experiment. The question is: How is it possible that a similar derivation and prescription 
adopted for QCD does not lead to the new chiral invariant vacuum advocated in the 
SUSY caseS? 

Our view of this problem is that we believe that in both cases (supersymmetric and 
non-supersymmetric) this suspicious state does not appear as a well-defined vacuum state; 
however in the SUSY case the situation might be less clear than in the ordinary QCD, see 
below. The argument is the following. Formally, such a state seemingly does appear as a 
solution of the equation of motion stemming from the effective Lagrangian. However, it 
does not appear to be a state stable against quantum fluctuations. Indeed, the modulus 
of the "order parameter" = (AA) vanishes in this candidate "vacuum". However, the 
phase of <fi is ill-defined at this point, and the matrix of second derivatives describing 
quantum fluctuations is also ill-defined. An additional source of ambiguity of the matrix 
of second derivatives is a freedom of redefinitions \<j>\ — > \<p\ n or even \<p\ — > log(|0|) in the 
effective Lagrangian, which are able to change the sign of Wfh for the candidate "vacuum" 
|0| = 0. This freedom of redefinitions of \<p\ is due to our lack of knowledge of the kinetic 
term in the effective SUSY Lagrangian, which otherwise would fix the correct definition 
of canonical field. Such a behavior of the effective potential is an indication that this state 
is not a genuine vacuum state stable against quantum fluctuations. In fact, the situation 
for the SUSY theories is a little bit more controversial because one could argue that in 
the SUSY case the vacuum energy is zero and thus all quantum fluctuations should cancel 
out, no matter what the eigenvalues for specific fluctuations are. Nevertheless, we believe 
that an accurate analysis of quantum fluctuations in the SUSY case in the background of 
the chiral invariant vacuum state within some suitable regularization would demonstrate 
that this state should be dismissed as an inappropriate candidate for the vacuum state. 

To end up this section, we would like to discuss an apparent problem related to the 
effective potential (|24l) . Proceeding by analogy with SUSY theories, we might naively 

8 We thank the Referee for bringing our attention to this potential problem. 
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expect that the effective potential, obtained for large m — > oo, could be continued to 
small masses m < A, where the heavy "glueball" fields h, h could be integrated out. (This 
is e.g. how the Affieck-Dine-Seiberg superpotential [Q was obtained in Ref.Q] from the 
VY effective potential for SUSY gluodynamics.) However, such a procedure gives correct 
results for supersymmetric theories just due to specific Ward identities which allow one to 
prove that the dependence of the gluino condensate on the chiral superfield mass m, viz. 
(AA) ~ ■sfm for SU(2) gauge group, obtained for SQCD for small m <C A, is actually exact 
and valid also for large m — > oo [1C[. As no such relation holds in the non-supersymmetric 
case, we actually have no reason to continue the above formulas to the region of small 
m ~ A. If we still do so, one can easily see that we do not reproduce in this way the 



anomalous term in the effective chiral Lagrangian of Ref . [ 14 1 . This seeming problem with 



Eq.([24]) is resolved by the fact that the "QCD limit" m -> and the "YM limit" m -> oo 
are actually separated by a kind of phase transition which changes the number of vacua 
in the theory. As the number of vacua is determined by the /3-function (see below in Sect. 
4), it should not change when a very heavy fermion is added to pure YM theory, and the 
effective theory below the decoupling scale Mq ~ m is considered. (This requirement will 
be given a formal content in the next section.) On the contrary, when the light quarks 
are present, the /3-function changes, and so does the number of vacua. Thus, analytic 
continuation of the above formulas to the small m < A region leads to a theory other 
than the one described by the effective chiral Lagrangian of Ref. ||14||. As is shown in 
|T5| , the latter is correctly reproduced by a different procedure. Namely, one should start 
directly from an effective Lagrangian for QCD with light quarks, which realizes at the 
tree level the anomalous conformal and chiral symmetries of QCD. The "glueball" part 
of this effective Lagrangian for QCD |15j is similar to Eg-fllCf). 



4 Global quantization, holomorphy, and Ward iden- 
tities in the effective Lagrangian approach 

In this section we would like to discuss a few related topics. First, we analyse the construc- 
tion of the effective potential ( p4[ ) of the u-theory while keeping track of multi-valuedness 
of effective potentials for both d- and u-theories. As we will see shortly, this analysis is 
consistent with the result (|28|) . Another purpose of this section is to discuss holomorphy 
of the effective theory in the fermion mass m from the point of view of the underlying 
QCD Lagrangian. This will allow us to establish the correspondence with an alternate 
method which was suggested by us earlier in Ref.|]lTJ to find the parameter £, where the 



same value £ = 4/(36) was obtained. It will be shown how the results of Ref. [11] follow 
from the effective Lagrangian approach. 

The analysis of the previous section was based on dealing with a fixed branch of 
the effective potential. As was discussed for SUSY theories in ||, and for pure YM 
theory in ||, the multi-valuedness of the effective potential necessitates the summation 
over all branches of a multi-valued action in the path integral (this is the prescription 
for constructing an improved effective potential). This procedure enforces some global 
quantization rules in the effective theory ||, [8| which reflect quantization of the topological 
charge in the fundamental theory. Technically, the global quantization in the d-theory 
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arises due to the Poisson formula 



(31) 

Let us now consider the problem of matching the d- and u- theories from the point of 
view of the integrating out procedure in the u-theory with account for the fact of multi- 
valuedness of the resulting potential. Given any branch of a multi-valued action, which is 
compatible with the Ward identities of the d-theory and the renormalization group, one 
achieves the single- valuedness of the partition function for the u-theory by the summation 
over all branches of the action in the path integral. In other words, we impose a universal 
rule for defining the partition function (improved effective potential) for both the d- and 
u-theories. Furthermore, it is natural to require that, after integrating out the field U in 
such a partition function for the u-theory, we should come back to the correct partition 
function of the d-theory. This requirement has a non-trivial content, since generically the 



logarithms in Eqs. (|12|) and (gj) imply different global structures in the complex h plane, 
i.e. different global quantization rules for the d- and u-theories. In this way we extend 
the integrating in/out procedure to match not only the perturbative scale redefinition in 
the low energy d-theory, but also the global, topological, properties of the latter. The 
meaning of this requirement is that turning on the heavy quark should not change 
the number of 9 vacuum sectors. As we will see shortly, this "topological" matching 
of the two theories agrees with the result ( |28| ) for the parameter £. 

Let us define the improved effective potential F for the u-theory in the same way as 
was done in Eq. fllQ|) . Using Eq. (|24]) , we obtain 



e 



-iVF 



Y< ex p\- iV 



n,m=—oc k,l 



l H 8 _ h 2 , _ 36 U 

4 {1 -Tb )hLo9 ^ + 3b hLo9 T^f 



2 q ( v 8 v 8 

+777 (h + h)-mU -mU-7r-lk-(l-—) + l- — + +n + m 
36 p \ q 36 q 36 




where the sums over k, I are finite. It is clear from Eq . (|32|) that the minimization condition 
F{j = leads to the same Eq. fl2"o]) . The reason for that is the absence of global quantization 
rules for the field U which is thus unconstrained. It can be readily seen that substituting 
this solution back to Eq. (|52D , we reproduce exactly the "dynamical" term in Eq. ([TDD • 
However, as we stated above, the "topological" term for the u-theory should also match 
the one of the d-theory. Let us now compare these two expressions. For the u-theory, the 
"topological" term is 

rr q / p . 8 , 7 p 8 9 \ h-h 

iV- K l[k^(l--) + l^-+n + m+ — ) — — , 33 
p \ q 36 q 36 lit J 1% 

where the integers n,m are ranged from (— oo) to (+oo), while the integers k, I reside on 
finite intervals. On the other hand, the "topological" term of the d-theory is 

iVn- Ik'- +n' + ^- ) ^--^ , n' = 0,±l,... , fc' = 0,l,...,g-l. (34) 
p \ q 27i J 2i 
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One can see that expression (|33|) has the same form as (151) only if 



- — = r — integer , (35) 
q 36 



for which expression (|33|) becomes 



8/36 „ , N 9 \ h-h 

> V7r Wr[ k j +n+ m + r{l ~ k) + s Hr ■ (36) 



Therefore, unless the constraint (|35| ) is satisfied, integrating out the field U in the u- 
theory does not lead back to the correct d-theory, but instead yields a theory with a 
different global quantization rule for the field h (i.e. different quantum mechanically). 
It is remarkable that these "topological" arguments are consistent with the value of £ 
which was arrived upon in Eq.(^) (still, they are less restrictive, as it usually happens 
for considerations based on topological arguments, than Eq.(p8[) above). On the other 
hand, in contrast to the argument of the previous section, no reference to the underlying 
theory was made above; only the internal consistency of the method just developed was 
used. The agreement between these two different lines of reasoning thus supports both 
the reliability of the calculated value of £, and self-consistency of the effective Lagrangian 
constructed. 

The next topic we would like to discuss is holomorphy of the effective Lagrangians. We 
have seen in the previous section that the holomorphic structure of the effective potential 
for the d-theory implies holomorphy in the fermion mass m of the effective potential for 
the u-theory. This means that m should be effectively treated as a complex parameter, 
with the correspondence rule ( p9|) between the operators of the effective and fundamental 
theories. Here we would like to discuss what such kind of holomorphy implies in terms of 
the fundamental QCD Lagrangian. 

Let us first note that in our case holomorphy refers to a physical heavy fermion with 
mass m — ■> oo. In the fundamental theory, the introduction of a heavy fermion of mass m 
requires a regularization on yet much higher ultra-violet scale Mr, which can be thought 
of as the mass of a Pauli-Villars regulator. In the infinite mass limit m, Mr — > oo the 
properties of the theory in respect to the physical and regulator fermions are identical 
(up to some sign differences) by definition. Thus, we may expect the fundamental theory 
to be holomorphic, in a sense, in the Pauli-Villars regulator mass Mr. The analysis of 
this paper shows what kind of holomorphy we may expect: it should be holomorphy 
of nonperturbative vacuum condensates or zero momentum correlation functions with 
perturbative tails subtracted, because these are the objects generated by the effective 
Lagrangian^. Precisely this kind of holomorphy for QCD with light quarks in the Pauli- 



Villars fermion mass Mr was suggested some time ago by Kuhn and Zakharov [IB| in 
a somewhat different context. Working in the chiral limit, these authors have related, 
using analyticity in Mr, the proton matrix element of the topological density (p\GG\p) 
to the matrix element (p\G 2 \p) which is fixed by the conformal anomaly. Recently, we 

9 As was mentioned in the Introduction, we do not expect that kinetic terms in the effective Lagrangian 
possess analogous holomorphic properties. Similarly, holomorphy is apparently of no use for correlation 
functions with non-zero momentum in the full theory. 
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have applied [JTTJ a similar idea to relate the zero momentum two-point function of GG 
to that of G 2 . The statement of holomorphy in the fermion mass m —>■ oo was rephrased 
there as a method to evaluate zero momentum correlation functions of the chiral fermion 
bilinears ^l^r and ^r^l with the operators G 2 ,GG from the known Ward identities 
involving the gluon operators only. As will be clear below, the value of £, the parameter of 
interest, follows from these relations. It is our purpose here to show that these results of 
Ref. [|y]] follow from the effective Lagrangian constructed in Sect. 3. Therefore, it comes 



as no surprise that in the effective Lagrangian approach we end up with precisely the 



same value of £ as that obtained in [fLl ]. (Again, we recall that the same regularization 
scheme is implied in the approach of [TTJ and the present paper.) 

Performing the inverse Legendre transform of the effective potential (1241) with respect 



to H and U ', we obtain the generating functional W m (J, j) of connected Green functions: 

1 /m+7'\ 8/(3b) 
W m (J,j) = — A cA* QCD l-^) e~* J . (37) 
4 \ A QCD j 

It is clear from the above analysis that the sources j, J are holomorphic: the differentiation 
in respect to j produces insertions of the chiral operator —^l^r, while the derivative 
d/dJ produces insertions of the operator H. As is seen from Eq. fl3~7|) , the differentiation 
with respect to j can be substituted by the differentiation with respect to m. In this way 
we find the correlation function 

% J dx(0\T{H(x) m* L * R (p)}\0) = -m^- ^ W m (J, 0) = ~(H) . (38) 

Analogously we can find further correlation functions: 

% j dx(0\T{H(x) mf A(0)}|0) = ~^(H) , (39) 

i J dx(0\T{H(x) m*i* B (0)}|0> = , (40) 
i J dx(0\T{H(x) mf r9 l (0)}\0) = . (41) 



The relations proposed in Ref. JTTJ now follow if we take linear combinations of Eqs. 
4"ip . Taking the sum of Eqs. (^3) and (^0|) and using 



H + H = -ba s /(8n)G 2 , (H) = (3b/2){m* L * R ) , 

we obtain 

i f dx (o\ T {^G 2 (x) m* L * fl (0)}|0) = ^■(mM R ) , (42) 

J 127T 60 

and analogously from Eqs. (|39|),(|4"lD 

i f rfx ( |T{-^G 2 (a;) m*A(0)}|0> = %r(mV R * L ) . (43) 

J lZ7T 60 

On the other hand, taking the differences for the same pairs of Eqs. (|38| - ^ID , we obtain 
i I dx(0\T{^GG(x) mV L * R (p)}\0) = £(mV L *R) , (44) 

J o7T 60 
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a. 



dx(0\T{-^GG(x) m* B *x(0)}|0) 

07T 



-i—{ m m R m L ) . 



(45) 



Eqs. ([42]-|4"5|) are precisely the relations suggested earlier in Ref.JTT| with a different 
motivation to find the topological susceptibility in terms of the gluon condensate. Taking 
the difference of Eqs. (fPj) and fl4"5]) and using Eqs.(p7|), we obtain 



(46) 



which is exactly the result obtained in JTT|. (As in [jTlJ], this scheme dependent result 



implies a particular regularization prescription where correlation functions are defined via 
the path integral with perturbative contributions subtracted, see the discussion above.) 
We have thus closed the circle: starting with Eq.(|3|) with an unspecified parameter £, we 
have found the value of £ using the effective Lagrangian, obtained from the anomalous 
Ward identities (|I]-fD, along with the integrating in/out technique, and eventually fixed 
the initial Eq.@. The introduction of a heavy quark was crucial to find the parameter £. 
One should note that the main assumptions of the present paper and that of Ref.[ll] are 



quite different, and yet lead to the same result ([46]). We consider this as an evidence in 
favor of correctness and self-consistency of our approach. 



5 Discussion and comparison with related works 

The purpose of this section is a qualitative discussion of the results obtained in the present 
paper. We will first address the counting of vacua following with our methods and compare 



its N c dependence with the results obtained for softly broken SUSY models |17j and the 
behavior found for the lattice Z p gauge models fl8| . Another purpose of this section is a 



qualitative comparison of our picture of the vacuum structure and 9 dependence in YM 



theory with a recent work by Witten [12] who studies the same issues in the limit N c — ► oo 



on the basis of the AdS / CFT correspondence [|T^] . 

As was discussed in ||, the value of £ determines the number of different nondegenerate 
9 vacuum sectors in pure YM theory. When the rational number p/q = l/(2£) is fixed 
with p and q being relatively prime, the number of different non-degenerate 9 vacua[^] is 
p. More precisely: 

number of different 9 sectors — p — 11 • N c for N c = odd 

A1N C 11N C UN Cl „ 
p = min.integery ; — - — ; J for N c = even , (47) 



where we have used £ = 4/(36). In particular, it follows from Eq.(47) that the number of 
vacua is 11 for N c = 2, and 33 for N c = 3. Although these numbers may look strange, 
they seem to be the only ones compatible with the effective Lagrangian and integrating 
in procedure considered in this paper. (They would be a wrong answer if one of the 
assumptions of the present approach were incorrect.) As was pointed out in [R], this 



10 Here we would like to recall that, for any generic value of 9, there is only one true physical vacuum 
which is a lowest energy state among all 9 vacua. For a fixed value of 9, additional vacua show up in the 
thermodynamics limit only indirectly via the value of parameter £ in the correlation function . 
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counting of vacua disagrees with what could be expected starting from SUSY YM theory 
broken softly by a small gluino mass m g <C A ||17|| . The latter theory predicts iV c vacua 
for small m g . Here is how it comes about (see a discussion by Shifman in [fj] for more 
detail). In the limit of small m g the VEV of the holomorphic combination G 2 + iGG is 
proportional to the VEV m g (XX) where the gluino condensate (AA) is to be calculated 
in the supersymmetric limit m g = 0. The 9 dependence of the latter is known flQfl : 



(AA) ~ exp(i9/N c + 2irk/N c ) , = 0, 1, . . . , N c — 1, which corresponds to N c degenerate 
vacua. When m g ^ 0, the vacuum degeneracy is lifted. For N c = 3 and 9 = 0, we have 
one state with negative energy E = — m 9 A| yM , and two degenerate states with positive 
energy E = (l/2)m 9 A| yM . The former is the true vacuum state of softly broken SUSY 
gluodynamics, while the latter are metastable states with broken CP. The lifetime of 
the metastable states is very large for small m g , and decreases as m g approaches A$ym 
. When 9 is varied, the three states intertwine, thus restoring the physical 2tt periodicity 
in 9. This picture implies the values p = N c , q — 1, different from those suggested by the 
present approach. 

The problem with the above SUSY-motivated scenario is that the genuine case of pure 
YM theory corresponds to the limit m g ^> A$ym which is not controlled in this approach. 
Thus, although naively one could expect that increasing of m g to higher values m g > A 
does not change the number of vacua of the theory, this expectation is unwarranted. It is 
conceivable that an additional level splitting occurs with passing the region m g ~ A SY m 
where the SUSY methods become inapplicable. Our results imply that this is indeed what 
happens, i.e. that there exists a sort of phase transition that separates the softly broken 
SUSY YM theory with m g <C A from pure non-supersymmetric gluodynamics. This 
expectation is consistent with numerous evidences within the soft SUSY breaking theories 
indicating that the naive decoupling limit m g ^> Asym produces results incompatible with 
the known infrared features of QCD. In particular, they include the wrong N c dependence 
of the effective chiral dynamics |2(| and the run-away behavior for Nf < N c in the presence 
of supersymmetry breaking pi]] . O n the other hand, it is curious to mention that the 
same value p = 11N C follows within a non-standard non-soft SUSY breaking suggested 
recently as a toy model to match the conformal anomaly of non-supersymmetric YM 
theory at the effective Lagrangian level. 

Our next remark concerns with another feature of the N c dependence in Eq. (fTFD . 
Naively, one could expect (and supersymmetric models support this expectation) that 
small variations of N c lead to small variations in the number of vacuum states. Our 
results suggest quite a different picture for certain values of N c : when the number of 
colors changes from N c = 8k — 1 to N c = 8k, the number of vacuum states abruptly 
changes from 11 (8A; — 1) to 11A;. At the same time, for generic values of N c , the variation 
of a number of vacua is quite smooth. Have we ever met with such kind of behavior in 
physics^]? The answer is yes. When a theory possesses two (or more) relevant integer 
parameters, the vacuum structure of the theory may undergo very dramatic changes with 
variations of these parameters. As an example one could consider Z p lattice gauge models 
in 4 dimensions in the presence of a 9 term which takes the values 9/(2n) = l/q where I 
and q are integers . In this case the physics is very sensitive to the numerical value of 



ii 



We thank the Referee for asking this question. 
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q, namely whether it is proportional to I or not. The physics changes drastically exactly 
at the points when q = I ■ n where n is another integer, much like in our case described 
above. (The concrete examples when 9/ (2tt) = 1/N and 9/(2it) = 2/N with N odd, were 
considered in detail in Ref. flSfl .) The physical explanation for such a behavior in the 
Z p model is based on the existence of a dual description where the relevant degrees of 
freedom are quite different from the original ones. In particular, exactly at these points 
one can construct a unique composite field which may condense according to different 
patterns depending on the value of 9/(2n). 

In a sense, the supersymmetric models are similar to the Z p lattice gauge models 
with 0/(2-71") — 1/N where a very smooth behavior is expected. We believe that the non- 
super symmetric models are closer to the case when 9/(2ir) = ^ in the Z p lattice gauge 
model with I, q being relatively prime. We do not know whether the explanation using 
the dual description for the Z p models can be extended to the case of gluodynamics, but 



an analogy with the examples discussed in Ref.|TB[ suggest that this might be the case. 

Finally, we would like to comment on another related development. Very recently, Wit- 
ten ]12| has shown how the qualitative features of the 9 dependence in non-supersymmetric 
YM theory - such as a multiplicity of vacua ~ iV c , existence of domain walls and exact 
vacuum doubling at some special values of 9 - can be understood using the AdS/CFT 
duality. The latter |19|] provides a continuum version of the strong coupling limit, with 
a fixed ultraviolet cutoff, for YM theory with N c —>■ oo, g\ M N c — > oo. As was shown in 
F2|, in this regime the 9 dependence of the vacuum energy in YM theory takes the form 



E vac (9) = C mm(9 + 2nk) 2 + 0(1/ N c 

k 



(4£ 



where C is some constant. We would like to make two comments on a comparison of our 
results with the picture advocated by Witten in the large N c limit. First, we note that 
the structure of Eq. ( ]4"8| ) agrees with our modified definition of the path integral including 
summation over all branches of a multi- valued (effective) action. Indeed, Eq. (|48|) suggests 
the correspondence 



C min [9 + 2nk) 2 lim log 

k V— >oo \ y 



e -vc(e+2wk) 2 

L k 



(49) 



using the definition of the vacuum energy through the thermodynamic limit of the path 
integral. With this definition which prescribes the way the volume V appears in the for- 
mula for the vacuum energy, the correspondence (^) appears to be the only possible one. 
On the other hand, the latter expression is exactly what arises (in the large N c limit) 
with our definition of the improved effective potential fllPl). Therefore, our prescription 
of summation over all branches of a multi- valued (effective) action seems to be consistent 
with the picture developed by Witten using an approach based on the AdS/CFT corre- 
spondence. In particular, our picture of bubbles of metastable vacua bounded by domain 
walls considered in the context of QCD within an effective Lagrangian approach in E7[ is 



in qualitative agreement with that suggested by Witten |TI| for the pure YM case. 

Second, one may wonder whether the approach of Ref. (T^j can provide an alternative 
way to fix the paramaters p, q of interest. We note that Eq. (|48|) indicates a non-analyticity 
at the values 9 C = ir (mod2^) only, where CP is broken spontaneously. If the technique 
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based on the AdS/CFT duality could be smoothly continued to the weak coupling regime 
of non-supersymmetric YM theory, this would result in the values q — 1, p ~ N c . However, 
the possibility of such extrapolation is unclear, as for small A = g\ M N c the background 
geometry develops a singular behavior and the supergravity approach breaks down. There 
might well be a phase transition |23[ when the effective YM coupling $y M N c is reduced. 
That such a phase transition should occur in the supergravity approach to QCD% was 
argued in |24| . Other reservations about the use of the supergravity approach to the non- 
supersymmetric YM theory in D=4 have been expressed in [25] where no perturbative 
indication was found for decoupling of unwanted massive Kaluza-Klein states of string 
theory. On the other hand, there exist some evidences from lattice simulations that a 
critical value of 9 moves from 9 C = it in the strong coupling regime to 9 C < it in the 
weak coupling regime []26[| . In terms of parameters p, q, such a case corresponds to q ^ 1. 
Therefore, we conclude that if no phase transition existed in the supergravity approach, 
our results would be in conflict with the latter which would imply p = 0(N c ),q = 1. 
In this case, the assumptions made in the present work would have to be reconsidered. 
Alternatively, there might be no conflict between the two approaches if such a phase 
transition does occur. 



6 Conclusions 

In this paper we suggested using the integrating in/out procedure to study the properties 
of a low energy Lagrangian for gluodynamics obtained in ||. We have shown that a 
particular holomorphic structure of this effective Lagrangian naturally corresponds to 
holomorphy in the fermion mass in an extended theory with a heavy fermion, provided the 
standard form of the fermion mass term is used. This observation supports the proposals 
of Refs. |I| and [[□]] where the idea of holomorphy in the regulator or heavy fermion 



mass was applied, respectively, to the study of matrix elements and correlation functions 
of the topological density operator GG. We have argued that the integrating in/out 
method provides not only a perturbative matching of the dimensional scale parameters 
in different theories, but also has a global, "topological", content. This holomorphic and 
"topological" matching of the theories with and without a heavy fermion was shown to 
fix, under certain assumptions, the number of different non-degenerate 9 vacua for YM 
theory. The result obtained implies that the number of vacua for softly broken SUSY YM 
gluodynamics changes discontinuously when the gluino mass becomes large. A check of 
this conclusion by different methods and analysis of its possible consequences would be 
an interesting problem for a further study. In particular, we have seen that the analysis of 
the 9 dependence is relevant for the question of presence or absence of a phase transition 
in the supergravity approach to the non-supersymmetric YM theory. In a more general 
context, it is perhaps worthwhile to point out that the modified path integral prescription 
of summation over all branches of a multi-valued (effective) action suggested for the 
case of pure YM theory in || (and previously proposed in different settings in Refs. 
|p8|,P) is consistent with the picture developed by Witten using an approach based on 
the AdS / CFT correspondence. Analogous modifications of the partition function for the 
case of lattice regularized sigma models and Abelian lattice models are found to be the 
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only self-consistent method [25] for the analysis of dualities on the lattice. This may 
provide further indications that such modifications of the partition function is the correct 
way to work with multi-valued actions. 



19 



References 

[1] N. Seiberg, Phys. Lett. B318 (1993) 469; Phys. Rev. D49 (1994) 6857. 

[2] K. Intriligator and N. Seiberg, |hep-th/9509~09~6"l 
M. Peskin, |hep-th/ 9 702094] 



M. Shifman, [hep-th/ 9704114 . 

[3] G. Veneziano and S. Yankielowicz, Phys. Lett. 113B (1982) 231. 

T. Taylor, G. Veneziano and S. Yankielowicz, Nucl. Phys. B218 (1983) 439. 

[4] K. Intriligator, R.G. Leigh, and N. Seiberg, Phys. Rev. D50 (1994) 1092. 
K. Intriligator, Phys. Lett. B336 (1994) 409. 

[5] I. Affleck, M.Dine, and N. Seiberg, Nucl. Phys. B241 (1984) 493; Nucl. Phys. B256 
(1985) 557. 

[6] L. Intriligator and N. Seiberg, Nucl. Phys. B431 (1994) 484. 
[7] D. Finnell and P. Pouliot, Nucl. Phys. B453 (1995) 225. 

I. Halperin and A. Zhitnitsky, Phys. Rev. D58 (1998) in press; |hep-ph/ 97 11398 . 



[9] A. Kovner and M. Shifman, Phys. Rev. D56 (1997) 2396; |hep-th/9702T74 . 
[10] M.A. Shifman and A.I. Vainshtein, Nucl. Phys. B296 (1988) 445. 
[11] I. Halperin and A. Zhitnitsky, Mod. Phys. Lett. A13 (1998) 1955; [hep-ph/ 9707286 . 
[12] E. Witten, |hep-th/9807109 . 



[13] V.A. Novikov, M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B191 
(1981) 301. 

[14] E. Witten, Ann. Phys. 128 (1980) 363. 

P. Di Vecchia and G. Veneziano, Nucl. Phys. B171 (1980) 253. 

[15] I. Halperin and A. Zhitnitsky, |hep-ph/9803301| . 

[16] J.H. Kiihn and V.I. Zakharov, Phys. Lett. B252 (1990) 615. 

[17] A. Masiero and G. Veneziano, Nucl. Phys. B249 (1985) 593. 

[18] J. Cardy and E. Rabinovici, Nucl. Phys. B205 (1982)1. 
J. Cardy, Nucl. Phys. B205 (1982)17. 

[19] J. Maldacena, |hep-th/9711200] 



S. Gubser, I.R. Klebanov and A.M. Polyakov, |iep-th/9802T09 
E. Witten, |hep-th/ 98021501 |hep-th/ 98031311 . 



[20] S. Martin and J. Wells, |hep-th/9801157 . 



20 



[21] J.L.F. Barbon and A. Pasquinucci, |hep-th/ 9804029" 



[22] F. Sannino and J. Schechter, |hep-th/ 9708111 . 
[23] D. Gross and H. Ooguri, |hep-th/9805I29". 



[24] J. Greensite and P. Olesen, |hep-th/ 9806235 



[25] H. Ooguri, H. Robins and J. Tannenhauser, [hep-th/9806171 



[26] G. Schierholz, Nucl. Phys. Proc. Suppl. 42 (1995) 270; |riep-lat/9412083| . 

A.S. Hassan, M. Imachi, N. Tsuzuki and H. Yoneyama, Prog. Theor. Phys. 95 (1996) 
175. 

H. Yoneyama, talk at Lattice-98. 
[27] T. Fugleberg, I. Halperin and A. Zhitnitsky, to appear. 
[28] A.V. Smilga, Phys. Rev. D49 (1994) 6836. 
[29] S. Jaimungal, |hep-th/98052lTi |iep-th/98080T8 . 



21 



